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ABSTRACT 


- —  y Four  regional  indices  of  sea  ice  variability  and  thirty  years  of  monthly 
sea  surface  temperature  (SST)  and  sea  level  pressure  (SLP)  data  are  used  to 
evaluate  the  large-scale  interactions  within  the  ocean-ice-atmosphere  systems 
of  the  North  Pacific  and  the  North  Atlantic.  A  direct  association  between 
positive  (negative)  SST  anomalies  and  light  (heavy)  ice  in  the  Bering  Sea  is 
indicated.  The  SST/ice  coupling  in  the  northern  Alaskan  waters  is  also 
statistically  significant.  The  SST  anomalies  in  both  the  North  Pacific  and 
the  North  Atlantic  appear  as  persistent  responses  to  atmospheric  forcing. 

The  SLP  distributions  over  both  oceans  correlate  significantly  with  the 
arctic  SLP  distribution,  although  the  midlatitude  atmospheric  circulation 
shows  no  systematic  tendency  to  lead  or  lag  the  arctic  circulation.  Sea  ice 
fluctuations  correlate  most  highly  with  the  atmospheric  forcing  over  ante¬ 
cedent  periods  of  approximately  9  months  in  the  northern  Alaskan  waters,  3-4 
months  in  the  Bering  Sea,  and  5-7  months  in  the  East  Greenland  Sea.  The 
strength  of  the  9-month  mean  onshore  air  flow  component  describes  over  70% 
of  the  northern  Alaskan  ice  variance.  + - 
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1.  INTRODUCTION 

A  sizable  literature  on  the  subject  of  large-scale  interactions  between 
sea  ice  and  the  atmosphere  has  evolved  in  recent  years.  The  existence  of 
coherent  fluctuations  of  ice  extent  over  thousands  of  kilometers  has  been 
identified  by  Walsh  and  Johnson  (1979a),  Lemke  et  al.  (1980)  and  Smirnov 
(1980).  Associations  between  interannual  ice  variability  and  fluctuations 
in  meteorological  fields  (e.g.,  sea  level  pressure,  surface  temperature,  700 
mb  height)  have  been  evaluated  in  varying  degrees  of  detail  by  Lamb  and  Morth 
(1978),  Zakharov  and  Strokina  (1978),  Rogers  (1978),  Rogers  and  van  Loon 
(1979)  and  Walsh  and  Johnson  (1979b).  While  modeling  experiments  have  shown 
statistically  significant  atmospheric  responses  to  prescribed  extremes  of 
sea  ice  extent  (Herman  and  Johnson,  1979),  the  consensus  emerging  from  the 
data-based  studies  is  that  large-scale  ice  anomalies  correlate  more  highly  with 
antecedent  than  with  subsequent  meteorological  fluctuations.  The  percentage 
of  ice  variability  that  may  be  specified  in  a  hindcast  sense  from  meteorologi¬ 
cal  data  appears  to  be  on  the  order  of  10-70%  depending  on  the  region,  the 
temporal  resolution,  the  quality  of  the  ice  data  and  the  choice  of  the 
meteorological  variables. 

Commonly  cited  sources  of  the  "unexplained"  ice  variance  are  oceanic 
factors  such  as  variable  ocean  currents  and  sea  surface  temperature  (SST) 
distributions  (e.g.,  Haupt  and  Kant,  1976).  Schell  (1970)  and  Rogers  and 
van  Loon  (1979)  have  demonstrated  the  compatibility  of  large-scale  SST  and 
sea  ice  variations  in  selected  regions  and  selected  years.  Considerably 
more  effort  has  been  expended  recently  in  the  evaluation*  of  the  coupling 
between  large-scale  SST  and  atmospheric  circulation  anomalies  over  the  North 
Pacific  (e.g.,  Namias,  1976;  Davis,  1978;  Elsberry  and  Raney,  1978;  Broccoli 
and  Harnack,  1981)  and  the  North  Atlantic  (e.g.,  Ratcliffe,  1973;  Rowntree, 
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1976;  Haworth,  1978).  Since  the  results  of  these  studies  have  indicated  the 
existence  of  statistically  significant  associations  between  SST  anomalies  and 
the  atmospheric  circulation,  the  link  between  SST  fields  and  the  high-latitude 
sea  ice  distribution  must  be  viewed  in  the  context  of  Figure  1.  The  "direct 
association",  in  which  high  SSTs  enhance  melting  and  retard  freezing,  is  to 
be  distinguished  from  the  "indirect  association",  in  which  the  sea  ice 
distribution  is  influenced  by  SST  through  an  SST-induced  effect  on  the 
atmospheric  circulation.  The  evaluation  of  the  SST-ice  coupling  from 
observational  data  is  complicated  by  the  possibility  that  both  coupling 
mechanisms  may  act  simultaneously. 

The  present  work  is  an  attempt  to  identify  the  associations  between 
the  SST  distribution,  the  atmospheric  circulation  and  the  interannual 
variability  of  sea  ice  in  several  regions  within  the  marginal  ice  zone:  the 
northern  Alaskan  waters  from  Cape  Prince  of  Wales  to  Demarcation  Point 
(summer),  the  Bering  Sea  (winter),  the  East  Greenland  Sea  (summer)  and  the 
Denmark  Strait  (winter).  Alternatively,  by  using  the  distribution  of  sea 
level  pressure  (SLP)  to  represent  the  low-level  atmospheric  circulation,  we 
wish  to  evaluate  the  relative  strengths  and  directions  of  the  three  links 
represented  in  Figure  1:  SST/ice,  SST/SLP  and  SLP/ice,  Such  an  evaluation 
can  have  direct  implications  about  causal  mechanisms  in  the  three-component 
system.  For  example,  a  cross-correlation  function  that  is  asymmetric  about 
0  lag  can  aid  in  the  establishment  of  cause  and  effect.  Cross-correlation 
functions  will  therefore  be  evaluated  in  an  attempt  to  answer  the  following 
questions: 

(i)  Are  high-latitude  SST  fluctuations  generally  causes  or  effects 
of  large-scale  sea  ice  fluctuations? 
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(ii)  Are  middle-  and  high-latitude  SST  fluctuations  generally  causes 
or  effects  of  atmospheric  circulation  anomalies? 

(iii)  How  does  the  link  between  SST  and  sea  ice  compare  quantitatively 
to  the  link  between  SLP  and  sea  ice? 

.'ithin  the  context  of  the  three-component  analysis,  several  additional  items 
/ill  be  addressed.  Among  these  are  the  seasonality  inherent  in  (i)-(iii), 
ifferences  between  the  air-sea-ice  systems  of  the  North  Atlantic  and  the 
,orth  Pacific,  and  the  time  scales  that  are  pertinent  to  the  evolution  of 
arge-scale  ice  anomalies.  The  latter  item  will  have  direct  implications 
f  or  the  predictability  of  large-scale  sea  ice  fluctuations.  Finally,  the 

* 

ixamination  of  the  indirect  association  of  Figure  1  will  require  the  evalu- 

/  ' 
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;  ition  of  "teleconnections"  or  lead/lag  relationships  between  the  arctic  and 
;  !  ddlatitude  atmospheric  circulations.  While  the  existence  of  significant 

ead/lag  relationships  of  this  type  would  assign  directionality  to  arctic/ 
^lidlatitude  interactions,  the  subject  of  arctic/midlatitude  teleconnections 
}ias  received  surprisingly  little  attention  since  Namias'  (1958)  analysis 
■jf  several  examples  of  arctic  circulation  anomalies. 

The  data  analysis  is  based  on  monthly  grids  of  SLP,  SST  and  sea  ice 
^xtent.  The  use  of  the  monthly  time  unit  places  some  constraints  on  the 
inferences  concerning  the  time  scales  of  the  system.  However,  the  available 
,lata  on  SST  and  sea  ice  preclude  more  detailed  temporal  resolution.  Spatial 
/ariability  will  be  incorporated  into  the  analysis  through  the  use  of  the 
iorainant  patterns  of  SST  and  SLP  variability.  Pattern  analysis  is  used 
oecause  SST  and  SLP  fields  are  characterized  by  large  atnounts  of  spatial 
..oherence  (e.g.,  Namias  and  Born,  1970;  Davis,  1976).  Considerable  redun¬ 
dancy  would  therefore  result  from  the  use  of  individual  grid  point  values. 


Section  2  describes  the  sources,  resolution  and  some  limitations  of  the 
data.  The  distributions  and  patterns  that  form  the  basis  of  the  correlative 


analysis  are  presented  in  Section  3,  while  the  results  and  implications  of 
the  correlative  analysis  are  discussed  in  Section  4. 

2.  DATA 

The  study  is  based  on  monthly  SLP  and  SST  data  grids  for  the  years 
1948-1977  and  on  monthly  sea  ice  grids  for  the  years  1953-1977.  A  total  of 
300  grids  of  sea  ice  and  360  grids  each  of  SLP  and  SST  are  therefore  used  in 
the  analysis. 

The  sea  level  pressure  data,  which  are  in  the  form  of  5°x5°  latitude- 
longitude  grids  from  20°N-90°N,  were  obtained  from  the  National  Center  for 
Atmospheric  Research  (NCAR) .  While  the  NCAR  grids  are  based  on  several 
sources  (Jenne,  1975),  the  primary  source  during  the  1948-1977  period  was 
the  U.S.  Navy  Fleet  Numerical  Weather  Center.  The  use  of  the  high-latitude 
portions  of  these  grids  introduces  the  risk  of  data  unreliability  (e.g. , 
Albright,  1980).  Trenberth  and  Paolino  (1980)  have  found  that  discontinuities 
associated  with  erroneously  high  arctic  pressures  appeared  primarily  in  the 
1920’s,  after  which  time  there  was  little  evidence  of  such  errors.  In 
addition,  tabulations  of  arctic  drift  station  data  were  used  in  the  present 
study  to  check  the  accuracy  of  the  NCAR  analyses  in  the  Arctic.  The 
discrepancies  between  the  observed  (monthly)  station  pressures  and  the 
analyzed  values  were  generally  less  than  2  mb.  The  errors  are  therefore 
considerably  smaller  than  the  typical  standard  deviations  (3-6  mb)  of  the 
high-latitude  pressures. 

The  sea  ice  grids  are  the  monthly  ice  concentration  grids  described  by 
Walsh  and  Johnson  (1979a).  These  grids  contain  the  fractions  (tenths)  of 


ice  coverage  regardless  of  ice  type  and  thickness  in  a  set  of  grid  squares 
60  n  mi  (111  km)  on  a  side.  As  discussed  in  Section  3,  the  ice  variables 
used  in  the  present  work  are  departures  from  the  25-year  monthly  mean  areal 
coverage  of  ice  in  several  longitudinal  sectors  of  30°-40°  each.  Each  such 
departure  was  obtained  by  a  summation  of  the  products  of  the  ice  concentra¬ 
tions  and  corresponding  ocean  areas  in  all  grid  squares  of  the  longitudinal 
sector. 

The  sea  surface  temperature  grids  were  obtained  from  several  sources. 

The  North  Pacific  data  are  in  the  form  of  monthly  5°x5°  latitude-longitude 
grids  covering  the  ocean  areas  from  20°N-60oN  and  130°E-110°W.  These  grids 
were  compiled  at  the  Scripps  Institute  of  Oceanography  and  were  supplied  to 
us  by  R.  Born  and  J.  Namias.  While  the  Scripps  data  set  covers  the  period 
from  January,  1947  through  the  present,  only  the  grids  for  1948-1977  are  used 
here.  North  Atlantic  grids  for  the  years  1948-1967  were  obtained  from  a 
"quality-controlled"  version  of  the  Tape  Data  Family  11  (TDF11)  which  was 
compiled  from  several  sources  by  the  National  Weather  Records  Center, 
Asheville,  NC.  These  grids  cover  the  ocean  areas  of  the  northern  hemisphere 
between  0®  and  80°N  with  a  resolution  of  5°x5®.  The  quality  control 
process,  in  which  large  values  of  normalized  monthly  departures  were  dis¬ 
carded,  was  carried  out  at  the  British  Meteorological  Office  in  Bracknell, 
England.  A  copy  of  the  "polished"  data  set  was  obtained  from  NCAR.  North 
Atlantic  SST  data  for  the  years  1963-1977  were  obtained  from  the  hemispheric 
surface  temperature  analyses  of  the  U.S.  Navy.  These  analyses  were  interpo¬ 
lated  bilinearly  from  the  Navy’s  63x63  polar  stereographic  grid  onto  the 
latitude-longitude  grid  of  the  TDF11  data  set.  The  North  Atlantic  grids 
from  each  source  were  then  converted  to  departures  from  the  monthly  mean 
SSTs  derived  from  that  source  only.  The  subtraction  of  monthly  means 


computed  from  the  individual  data  sets  served  to  minimize  the  effect  of  any 
discontinuities  introduced  by  the  consolidation  of  the  two  data  sets. 
Correlations  between  the  two  sets  of  monthly  anomalies  were  computed  for 
the  five  years  (1963-1967)  of  data  overlap;  these  correlations  ranged  from 
0.70  to  0.92  over  a  set  of  50  grid  points  in  different  portions  of  the  North 
Atlantic.  The  common  variance  between  the  two  data  sets  therefore  ranged 
from  0.49  to  0.83.  Since  the  discrepancies  represent  substantial  portions 
of  the  SST  variance,  the  results  of  a  correlative  analysis  based  on  the 
North  Atlantic  SST  data  must  be  regarded  as  lower  bounds  on  the  coupling 
between  the  SST  fields  and  the  other  variables.  As  a  compromise  between  the 
uncertainties  in  the  two  data  sets,  algebraic  means  of  the  monthly  anomalies 
from  the  two  sources  were  used  for  the  1963-1967  period. 

It  should  be  noted  that  the  Scripps  and  TDF11  data  sets  often  do  not 
include  data,  values  within  5°  latitude  of  the  ice  edge.  While  these  missing 
values  are  attributable  to  the  scarcity  of  observations  near  the  ice  edge, 
they  limit  the  evaluation  of  the  "direct"  association  between  SST  and  ice 
extent.  The  analyses  in  the  following  two  sections  are  based  on  50  North 
Pacific  and  50  North  Atlantic  grid  points  that  were  chosen  for  spatial 
uniformity  subject  to  the  constraint  that  no  single  point  be  missing  from 
more  than  15%  of  the  total  number  of  grids. 

3.  DATA  REPRESENTATIONS 

a.  Sea  ice 

The  areal  coverage  of  sea  ice  was  evaluated  for  Atlantic  and  Pacific 
sectors  of  30-40°  longitude  for  three-month  periods  centered  on  the  approxi¬ 
mate  months  (March  and  August)  of  minimum  and  maximum  ice  extent.  The 
longitudinal  sectors  will  be  referred  to  as  the  northern  Alaskan  waters 
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(NA,  140°W-170°W) ,  the  Bering  Sea  (BS,  150°W-180°W) ,  the  East  Greenland  Sea 
(EG,  20°W-20°E)  and  the  Denmark  Strait  northwest  of  Iceland  (DS,  10oW-40°W) . 
Ice  extent  is  presented  as  sums  of  departures  from  the  normal  areal  coverage 
for  July-September  in  the  northern  Alaskan  and  East  Greenland  regions,  and 
as  the  corresponding  sums  for  February-April  in  the  Bering  Sea  and  the 
Denmark  Strait.  The  analysis  therefore  addresses  wintertime  ice  variability 
in  two  southern  regions  and  summertime  ice  variability  in  two  northern 
regions  within  the  seasonal  sea  ice  zone.  Figure  2  shows  the  25-year  time 
series  of  ice  coverage  in  each  of  the  four  regions.  Considerable  interannual 
variability  is  apparent  in  every  region,  although  the  two  plots  for  the 
winter  season  (IgS  and  I  )  show  a  greater  tendency  for  the  clustering  of 
heavy  and  light  ice  years.  Ice  extent  in  the  Bering  Sea  was  generally  below 
the  25-year  mean  during  the  1960's  and  above  the  mean  during  the  1970's, 
while  ice  in  the  Denmark  Strait  was  generally  below  normal  during  the  first 
half  of  the  period  and  above  normal  during  the  1965-1970  period.  However, 
none  of  the  25-year  trends  were  found  to  be  statistically  significant.  The 
increase  in  Bering  sea  ice  during  the  early  and  middle  1970 's  was  also 
reported  by  Niebauer  (1980) .  The  time  series  for  the  northern  Alaskan  ice 
correlates  at  -0.895  with  Barnett's  (1980)  Alaskan  ice  severity  index,  which 
ranges  from  0  for  the  most  severe  shipping  season  (1975)  to  624  for  the  most 
favorable  shipping  season  (1958)  of  the  25-year  period. 

Table  1  contains  the  correlations  between  the  four  regional  indices  of 
ice  coverage.  The  0-lag  correlations,  which  are  based  on  pairs  of  values 
from  the  same  calendar  year,  were  found  to  be  positive  ft>r  regions  in  the  same 
hemisphere  (Atlantic  or  Pacific)  and  negative  for  regions  in  opposite  hemi¬ 
spheres.  Since  the  95%  significance  level  for  a  sample  of  25  (independent) 
pairs  of  values  is  0.40,  the  EG/DS  correlation  of  0.57  indicates  a 


statistically  significant  tendency  for  winter  anomalies  in  the  Denmark  Strait 
to  be  followed  by  summer  anomalies  of  the  same  sign  in  the  East  Greenland 
Sea.  The  negative  values  of  the  Atlantic  vs.  Pacific  correlations  (EG  vs. 

NA,  EG  vs.  BS,  DS  vs.  BS)  provide  modest  quantitative  support  for  Smirnov's 
(1980)  contention  that  ice  anomalies  in  the  eastern  and  western  Arctic  are 
generally  of  opposite  sign.  These  negative  correlations  are  also  consistent 
with  recently  computed  sets  of  sea  ice  eigenvectors  (Lemke  et  al. ,  1980; 

Walsh  and  Johnson,  1979b)  which  indicate  that  the  most  common  modes  of  large- 
scale  ice  variability  contain  positive  anomalies  in  some  longitudes  and 
negative  anomalies  in  others.  The  lag  0  correlations  between  the  Atlantic 
and  Pacific  indices  are,  however,  only  marginally  significant  at  the  95% 
level. 

The  lag  1  correlations  of  Table  1  show  the  tendency  for  departures  in 
the  Bering  Sea  and  the  Denmark  Strait  to  be  of  the  same  sign  as  the  departures 
during  the  previous  summer  in  the  same  hemisphere.  In  addition,  the  year-to- 
year  persistences  (0.50  and  0.52)  of  the  two  wintertime  series  contrast  with 
the  negligible  persistences  (0.03  and  -0.07)  of  the  summer  series.  The 
relatively  large  persistences  of  the  wintertime  indices  are  statistical 
measures  of  the  tendency  for  the  clustering  of  heavy  and  light  ice  years  in 
the  two  wintertime  series  (Igg  and  ID£.)  of  Figure  2. 
b.  Sea  level  pressure  and  sea  surface  temperature 

In  view  of  the  large  volume  of  SLP  and  SST  data,  much  of  the  analysis 
and  discussion  is  based  on  the  dominant  patterns  of  variability  of  the  data 
fields.  These  patterns  are  empirical  orthogonal  functions  (EOFs),  which  are 
also  referred  to  as  eigenvectors  or  principal  components.  The  first  eigen¬ 
vector  of  a  set  of  data  fields  is  the  dominant  mode  of  variability  in  the 
sense  that  it  describes  a  greater  percentage  of  the  variance  of  the  data 


fields  than  does  any  other  pattern.  Each  succeeding  eigenvector  describes 
a  maximum  of  the  variance  that  is  unexplained  by  the  previous  eigenvectors. 

The  coefficient  or  amplitude  of  an  eigenvector  is  a  measure  of  the  extent  to 
which  that  eigenvector  pattern  is  present  in  a  particular  anomaly  field. 

The  various  time  series  of  the  eigenvector  coefficients  form  the  basis  of 
the  statistical  analysis  in  Section  4.  Among  the  advantages  of  the  eigen¬ 
vector  representation  are  the  effective  compression  of  the  data  and  the 
orthogonality  (independence)  of  the  patterns  in  space  and  time. 

Since  a  goal  of  the  analysis  was  the  evaluation  of  the  directions  and 
relative  strengths  of  the  links  in  Figure  1,  sets  of  eigenvectors  were 
computed  for  SST  in  the  midlatitude  ocean  and  SLP  in  middle  and  high  latitudes. 
These  three  sets  of  patterns  were  computed  separately  for  the  North  Pacific 
sector  and  for  the  North  Atlantic  sector.  All  eigenvectors  were  computed 
using  normalized  departures  from  monthly  means.  The  normalization  involved 
a  subtraction  of  the  monthly  mean  (SLP  or  SST)  for  the  grid  point  and  a 
division  by  the  standard  deviation  for  that  grid  point  and  calendar  month. 

This  procedure  removes  the  seasonal  cycle  that  otherwise  would  concentrate 
most  of  the  variance  into  a  single  mode  (e.g.,  Weare,  1977)  and  avoids  an 
overweighting  of  high-variance  regions  in  the  resulting  patterns.  Because 
the  spatial  domain  of  each  set  of  eigenvectors  in  the  present  study  is 
relatively  limited,  the  spatial  bias  introduced  into  the  eigenvectors  by  the 
normalization  is  less  problematic  than  in  studies  of  hemispheric  patterns 
(e.g.,  Rogers,  1979). 


Figures  3  and  4  show  the  first  four  eigenvectors  of  the  Pacific  SLP  sets, 

MP  HP  1 

P  (Midlatitude  pacific)  and  P  (High-latitude  Pacific) .  Each  set  was 

computed  from  a  network  of  50  grid  points  represented  by  dots  in  the  figures. 

Table  2  summarizes  the  major  features  of  the  various  sets  of  eigenvectors. 

MP  MP 

Both  P^  and  P^  contain  north-south  gradients  from  which  anomalous  east-west 

geostrophic  airflow  may  be  inferred.  The  strongest  pressure  gradient  is 

M  P  M  H 

approximately  5-10°  farther  southward  in  P^  than  in  P^  .  The  results  in 

MP  MP 

Section  4a  will  support  the  hypothesis  that  P^  and  P^  are  summer  and  winter 

manifestations  of  the  same  pattern,  which  corresponds  to  a  mass  exchange 

between  the  Aleutian  low  and  the  Pacific  subtropical  high.  Southward  shifts 

of  5-10°  from  summer  to  winter  are  typical  of  these  semi-permanent  features. 

MP  MP 

P^  and  P^  represent  somewhat  more  complex  anomalies,  although  the  implicit 

north-south  airflow  over  the  Bering  Sea  is  pertinent  to  the  interpretation 

KP  MP 

of  the  correlative  results  in  Section  4.  The  patterns  P,  -  P,  describe 

~1  ~4 

65.4%  of  the  variance  contained  in  the  entire  set  of  50  midlatitude  Pacific 
pressure  eigenvectors. 

HP  HP 

The  first  four  high-latitude  Pacific  pressure  patterns,  P^  -  P^  ,  are 

shown  in  Figure  4  and  summarized  in  Table  2.  These  four  patterns  describe 

HP 

78.6%  of  the  50-point  variance.  While  represents  a  general  excess  or 

HP  HP 

deficit  of  mass  with  little  gradient  over  the  entire  domain,  P ^  and  P^ 

describe  large  gradients  in  approximately  perpendicular  directions.  The 

HP 

corresponding  geostrophic  airflow  and  ice  drift  are  north-south  in  P  and 


1  MP  P 

In  the  sections  that  follow,  vector  notation  (P  ,  T  ,  etc.)  is  used  to 

MP  P 

denote  eigenvectors.  The  corresponding  scalar  symbols  (P  ,  T  ,  etc.) 


denote  time-dependent  amplitudes. 
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HP  HP 

east-west  in  P^  .  The  importance  of  P2  in  the  evolution  of  Alaskan  ice 

anomalies  will  become  apparent  in  Section  4. 

Comparisons  can  be  made  between  the  patterns  of  Figures  3-4  and  Rogers' 

(1979)  hemispheric  eigenvectors,  although  the  latter  were  computed  for  the 

MP  , 

winter  season  only.  P2  appears  to  correspond  with  Rogers'  second  unnormal¬ 
ized  eigenvector,  which  in  turn  resembled  the  so-called  "North  Pacific 

Oscillation"  in  pressure.  On  the  other  hand,  the  pattern  of  zonal  flow  in 
HP 

P^  is  quite  similar  to  the  high-latitude  circulation  in  Rogers'  first 
unnormalized  eigenvector,  which  showed  a  strong  resemblance  to  the  "North 
Atlantic  Oscillation"  or  "Seesaw"  (Rogers  and  van  Loon,  1979) .  Further 
comparisons  are  prevented  by  the  limited  spatial  domain  used  in  the  present 
work  and  by  the  restriction  of  the  Rogers  and  van  Loon  analysis  to  the  winter 


season. 

The  first  four  patterns  of  North  Pacific  sea  surface  temperature  are 

shown  in  Figure  5  (see  also  Table  2).  Since  these  patterns  describe  only 

48%  of  the  50-eigenvector  SST  variance,  we  may  infer  that  the  SST  fields 

contain  more  small-scale  variability  than  do  the  North  Pacific  SLP  fields. 

P  P  P 

Positive  amplitudes  of  T^  ,  T2  and  T^  correspond  to  negative  temperature 
-nomalies  in  the  Aleutian  region,  although  rather  large  positive  anomalies 

p 

are  found  south  of  45°N.  The  upper-air  anomaly  deduced  from  positive  T^ 
and  the  thermal  wind  relationship  is  an  enhanced  ridge  over  the  cential 
North  Pacific  with  an  anomalous  southward  component  of  geostrophic  airflow 
over  the  eastern  North  Pacific.  (Such  reasoning  is  based  on  the  assumption 
that  the  SST  anomalies  are  representative  of  the  temper&ture  gradients  in 
the  lower  troposphere  [Harnack  and  Broccoli,  1979]).  Similarly,  positive 


T2  implies  enhanced  zonal  flow  at  upper  levels  above  the  30°-40°N  latitude 

p 

belt.  The  relatively  strong  anomaly  center  in  T^  at  35°N,  160°W  is  very 
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close  to  the  center  (40°N,  170°W)  of  the  correlation  maps  used  by  Namias 

(1972)  in  his  analysis  of  the  spatial  scales  of  North  Pacific  SST  anomalies, 
p 

is  also  very  similar  to  Davis'  (1976)  first  eigenvector  of  North  Pacific 
SST. 

(2)  Atlantic  patterns 

MA  MA 

Figure  6  shows  the  first  four  patterns,  P,  -  P,  ,  of  North  Atlantic 

-L  4 

MA  MA 

sea  level  pressure.  P  and  P^  correspond  very  closely  to  the  first  two 

North  Pacific  patterns.  Variations  in  the  intensity  of  the  midlatitude  zonal 

flow  are  described  in  each  case.  The  third  North  Atlantic  pattern,  which  is 
...  MP 

very  similar  to  P  ,  represents  enhanced  northward  flow  in  the  North  Atlantic 
waters  east  of  Greenland.  The  first  four  North  Atlantic  patterns  describe 
60.2%  of  the  50-eigenvector  SLP  variance. 


The  first  four  patterns  of  high-latitude  North  Atlantic  pressure,  P^ 


HA 


„HA 


P^  ,  are  shown  in  Figure  7  and  summarized  in  Table  2.  A  close  correspondence 

.  _  '  HA 

with  the  Pacific  set  is  again  found  in  the  very  weak  gradient  of  P^  and  the 

north-south  and  east-west  gradients  of  P^  and  P^  ,  Although  the  double 

HA 

lobe  of  high  pressure  in  the  southern  portion  of  P„  is  somewhat  different 

J 

HP 

from  the  single  lobe  in  P^  ,  the  fourth  eigenvectors  of  the  two  sets  are 
extremely  similar.  The  first  four  high-latitude  North  Atlantic  patterns 
describe  76.9%  of  the  50-eigenvector  variance.  The  patterns  P?^  and  P^ 
are  both  quite  consistent  with  Rogers*  (1979)  first  hemispheric  eigenvector 
and  with  the  March  sea  level  pressure  difference  fields  associated  with  the 
North  Atlantic  Oscillation  or  Seesaw  (Rogers  and  van  Loon,  1979,  Figure  6). 
This  similarity  is  detectable  even  though  the  eigenvector  sets  P^  and  P^ 
are  computed  from  data  for  all  calendar  months. 

A  A 

Figure  8  contains  the  first  four  patterns  (T,  -  T. )  of  North  Atlantic 

-1  ">4 

A 

sea  surface  temperature,  T.  resembles  the  first  Pacific  SST  pattern  with  a 

f  X 


large  anomaly  in  the  central  portion  of  the  grid  and  weaker  anomalies  of  the 

A 

opposite  sign  in  the  northeast  and  southwest  corners.  If  is  representa¬ 
tive  of  temperature  anomalies  in  the  lower  troposphere,  then  the  thermal 
wind  constraint  implies  enhanced  ridging  above  the  lower  troposphere  over 
the  North  Atlantic.  ,  which  describes  anomalies  of  a  single  sign  over 

p 

all  but  the  west-central  portion  of  the  grid,  shows  some  resemblance  to  T^. 

The  generally  north-south  gradient  in  ,  on  the  other  hand,  is  similar  to 
P  A 

that  found  in  •  T^  contains  some  of  the  finer  structure  that  is  typical 

A  A 

of  the  higher  modes.  The  percentage  of  SST  variance  described  by  T^  -  T^, 

is  essentially  the  same  as  the  48%  described  by  the  first  four  North 
Pacific  SST  patterns.  These  four  patterns  correspond  quite  well  with  the 

high-latitude  portions  of  Weare's  (1977)  first  four  Atlantic  SST  patterns 

A  A  A 

if  the  order  of  T^  and  T^  is  reversed.  However,  the  agreement  between  T^ 

and  Weare's  second  eigenvector  is  somewhat  weak  in  the  eastern  North 

Atlantic. 


4.  CORRELATIVE  ANALYSIS 

The  results  of  the  analysis  will  be  described  separately  for  the  North 
Pacific  (Section  4a)  and  the  North  Atlantic  (Section  4b).  Pertinent 
similarities  and  differences  between  the  two  regions  will  be  noted. 
a.  North  Pacific 

Cross-correlations  between  the  amplitudes  of  the  dominant  patterns  of 
SST  and  middle-  and  high-latitude  SLP  were  computed  in  order  to  assess  the 
teleconnections  in  the  regional  ocean-atmosphere  system.  The  360  monthly 
grids  (1948-1977)  of  each  variable  were  used  in  the  computations.  The 
correlations  are  based  on  3-month  means  of  the  pattern  amplitudes  in  order 


to  remove  some  of  the  high-frequency  atmospheric  noise.  (While  the  magnitudes 
of  the  cross-correlaticns  were  found  to  increase  when  the  averaging  period 
was  made  longer  than  3  months,  this  increase  was  at  the  expense  of  temporal 
resolution).  A  two-dimensional  format  is  used  in  Figure  9  to  show  the 
dependence  of  several  of  the  correlations  on  season  (calendar  month)  and  on 
lag,  which  ranges  from  -6  to  +6  months  for  each  variable.  Stippling  is  used 
to  identify  those  portions  of  the  lag/season  domain  where  the  magnitudes  of 
the  correlations  exceed  0.4.  The  95%  and  99%  significance  levels  for  a 

sample  size  of  30  (years)  are  0.36  and  0.48. 

.  MP  P 

The  correlations  <P  T  >  between  North  Pacific  sea  level  pressure  and 

MP  P 

sea-surface  temperature  are  quite  large  for  several  pattern  pairs.  <~P^  T^> 
exceeds  0.8  during  the  winter,  while  the  magnitudes  of  several  other 
correlations  exceed  0.7  in  certain  portions  of  the  year.  Several  features 
of  the  <PMPT^- diagrams  have  important  implications  concerning  the  nature  of 
the  air /sea  coupling.  The  first  such  feature  is  the  asymmetry  about  0  lag 
in  the  distribution  of  high  correlations.  Since  a  large  majority  of  the 
correlations  exceeding  0.4  are  found  at  positive  lags  of  sea  surface  temper¬ 
ature,  the  coupling  is  evidently  one  in  which  the  atmosphere  tends  to  drive 

the  ocean.  Davis  (1976)  has  found  a  similar  directionality  in  his  North 

MP  P 

Pacific  analysis.  A  second  noteworthy  feature  of  the  <P  T  >  diagrams  is 
the  pronounced  tendency  for  the  shaded  regions  to  tilt  downward  and  to  the 
right.  This  tendency  suggests  that  the  directionality  noted  above  is  due 
primarily  to  the  persistence  of  SST  anomalies  after  an  oceanic  response  to 
atmospheric  (SLP)  forcing.  Elsberry  and  Raney  (1978)  have  shown  that  SST 
anomalies  at  individual  points  (Ocean  Weather  Ships)  are  attributable  largely 
to  limited  periods  of  strong  atmospheric  forcing  (e.g.,  storms).  The 
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Mr  r 

downward/rightward  tilt  of  the  larger  correlations  in  the  <P  T  >  diagrams 

of  Figure  9  may  be  manifestations  of  the  effects  of  this  type  of  atmospheric 

forcing  on  a  larger  spatial  scale  and  on  the  seasonal  time  scale.  Such 

interpretations  must  be  considered  speculative,  however,  until  supported  by 

more  detailed  analyses  of  the  evolution  of  SST  anomaly  fields. 

MP  P 

The  signs  of  the  <P  T  >  correlations  are  generally  consistent  with 

reasoning  based  on  advection  of  SSTs  by  SLP-induced  current  anomalies.  For 
MP  P  MP  P  MP  P 

example,  <P^  T^>  ,  <P^  T^>  and  <P^  T^>  represent  cases  where  negative 
(positive)  SST  anomalies  are  found  in  regions  of  anomalous  southward  (north¬ 
ward)  or  westward  (eastward)  transport  from  normally  colder  (warmer)  waters. 

MP  P 

The  interpretation  of  the  large  positive  <P^  T^>  is  less  straightforward. 

MP  P 

Since  the  coupling  between  P^  and  is  clearly  a  wintertime  phenomenon, 

MP  MP 

the  earlier  interpretation  of  P^  as  a  wintertime  counterpart  of  P^  is 

MP 

supported.  Positive  P^  corresponds  to  a  weakened  Aleutian  low,  while  a 

reduction  of  storm  activity  over  the  North  Pacific  is  consistent  with  a 

p 

shallow  mixed  layer  and  the  positive  SST  anomalies  of  T^.  Such  reasoning 

should  apply  primarily  in  the  late  winter  and  early  spring  (Elsberry  and 

Garwood,  1978),  although  Niebauer  (1980)  finds  a  significant  negative  . 

correlation  between  wintertime  cloudiness  and  annually  avevaged  SST. 

The  coupling  between  the  low-level  atmospheric  circulations  over  the 

HP  MP 

Arctic  and  the  North  Pacific  is  shown  by  the  <P  P‘  >  correlations  of 

Figure  9.  The  coupling  is  consistent  with  the  signs  of  the  anomalies  within 

HP  MP  HP  MP 

the  individual  eigenvectors,  e.g.,  positive  <P^  P ^  >,  negative  <P^  P  ^>, 

HP  MP  • 

positive  <P^  Pj  >•  However,  the  coupling  is  generally  weaker  than  the 

coupling  between  North  Pacific  SLP  and  SST  in  terms  of  both  the  magnitudes 


of  the  largest  correlations  and  the  maximum  lags  of  the  statistically 
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significant  correlations.  The  coupling  is  especially  weak  during  the  late 

HP  MP 

spring  and  early  summer.  A  pertinent  feature  of  the  <P  P  >  diagrams  is 
their  general  symmetry  about  0  lag.  There  is  therefore  no  apparent  tendency 
for  the  arctic  circulation  to  drive  or  be  driven  by  the  North  Pacific 
circulation. 

Finally,  the  coupling  between  the  arctic  SLP  and  North  Pacific  SST 

HP  P 

fields  is  _hown  by  the  <P  T  >  diagrams  of  Figure  9.  Not  surprisingly,  this 

coupling  is  even  weaker  than  the  Arctic/Pacific  SLP  coupling.  The  largest 

HP  P 

correlations  are  the  positive  wintertime  <P^  T2>*  imply  that  southward 

airflow  over  the  Bering  Strait  is  associated  with  below-normal  SSTs  over  a 

broad  area  of  the  North  Pacific  from  35°N  to  55°N.  The  correlations 

of  0.5-0. 6  agree  well  with  Niebauer's  (1980)  correlations  for  the 

HP  MP 

Pribiloff  Island  region  during  1973-1979.  As  in  the  case  of  the  <P  P  > 

correlations,  there  is  no  apparent  tendency  for  the  arctic  SLP  fields  to 

lead  or  lag  the  SST  fields  in  any  systematic  manner. 

The  relationships  between  summertime  ice  extent  along  the  northern 

P  HP  HP  MP 

Alaskan  coast  (I  Figure  2)  and  the  amplitudes  of  T  ,  P  ,  P_  and  P"^ 

WA  ~  j  *•  j 

are  shown  in  Figure  10a.  The  amplitudes  of  the  SLP  and  SST  patterns  were 

again  smoothed  with  a  3-month  running  mean  prior  to  the  computations.  The 

correlations  are  plotted  against  the  calendar  month  of  the  midpoint  of  each 

3-month  period.  Figure  10a  contains  plots  for  all  those  SLP  and  SST  variables 

that  correlate  with  I  at  a  statistically  significant  (95%)  level  in  at 

HP 

least  two  consecutive  3-month  periods.  (The  small  values  of  <IV,.P„  >  are 

NA  3 

HP 

plotted  only  for  comparison  with  >).  The  patterns  correlating  most 

WA  4 

P  HP  HP 

highly  with  1^  are  T^  and  Pj  .  A  relationship  between  1^  and  Is  not 

HP 

surprising,  since  positive  describes  a  flow  pattern  with  a  strong 
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onshore  component  in  the  northern  Alaskan  region.  However,  the  statistical 

HP 

significance  of  the  correlation  of  1  with  of  the  previous  winter  and 

spring  was  not  expected.  While  I  correlates  most  highly  with  the  May-June- 

HP  HP 

July  value  of  P2  (r  =  0.61),  1^  correlates  significantly  with  the  P2 

values  of  every  month  of  the  December-July  period  except  February.  The 

results  therefore  imply  that  onshore  airflow  during  the  previous  winter  and 

spring  does  contribute  to  the  severity  of  summer  ice  conditions  along  the 

northern  Alaskan  coast.  Later  figures  will  support  the  interpretation  of 

summer  ice  extent  as  an  integrated  response  to  atmospheric  forcing  over  the 

previous  6-12  months.  By  contrast,  the  alongshore  (east-west)  component  of 
HP 

the  airflow,  P^  ,  does  not  correlate  significantly  with  I  ^  during  any  month. 

p 

The  magnitudes  of  the  correlations  between  T-  and  T  .  are  also  unexpect- 

1  NA 

P  HP 

edly  large.  <Ij^T  >  exceeds  even  <INAP2  >  t*ie  3-month  periods  centered 

on  every  calendar  month  except  January.  The  high  persistence  of  SST  anoma- 

P 

lies  is  evident  in  the  slow  decay  of  <I„.T, >  after  the  summer  ice  season. 

NA  1 

p 

The  positive  sign  of  <I^T^>  implies  that  above-normal  ice  extent  along  the 

northern  Alaskan  coast  (I„.  >  0)  is  associated  with  above-normal  SSTs  in 

NA 

the  30°N-45°N  band  of  the  North  Pacific  and  with  below-normal  SSTs  in  the 
Aleutian  region  and  along  the  west  coast  of  North  America  (T^  >  0).  While 
the  mass  (water)  influx  through  the  Bering  Strait  is  northward  into  the 
Chukchi  Sea,  inflow  rates  (e.g..  Coachman  and  Aagaard,  1974)  appear  to  be 
too  low  to  permit  a  "direct"  association  between  winter/spring  Bering  Sea 
SST's  and  northern  Alaskan  ice  conditions  during  the  following  summer.  If 
Bering  Strait  inflow  is  a  factor  in  the  interannual  variability  of  northern 
Alaskan  ice  conditions,  current  fluctuations  and  attendant  variations  in 
water  stress  must  also  be  considered.  An  alternative  hypothesis  concerning 


A — 
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the  relatively  large  positive  values  of  cI^T  >  is  that  anomalously  strong 

southward  airflow  (positive  )  prior  to  heavy  ice  summers  is  a  factor 

contributing  to  both  heavy  Alaskan  ice  and  to  the  below-normal  SSTs  in  the 

Aleutian  region.  However,  this  interpretation  does  not  explain  the  fact 
P  HP 

that  <I^T  >  exceeds  <1  P2  >  nearly  every  month.  The  interplay  between 
surface  airflow,  ocean  temperatures  and  currents,  and  seasonally  varying  ice 
conditions  in  the  Alaskan  waters  is  undoubtedly  complex.  More  detailed  and 
systematic  measurements  of  oceanic  parameters  are  clearly  required  before 
the  interpretations  described  above  can  be  regarded  as  more  than  the 
speculative  rationalizations  of  observational  data. 

Figure  10b  shows  the  correlations  between  ice  extent  during  February- 

p  p  up  vrp 

April  in  the  Bering  Sea,  Ifis,  and  the  amplitudes  of  T^,  and 

HP 

P^  .  No  other  SST  or  SLP  fields  correlated  at  0.5  or  higher  with  1^  in  any 
month.  The  positive  correlations  between  l_c  and  the  first  two  SST  patterns 
indicate  an  association  between  heavy  Bering  Sea  ice  and  below-normal  SSTs 

p  p  p 

(T^  >  0  ,  T2  >  0)  in  the  Aleutian  region.  While  T2  correlates  significantly 

p 

with  XBS  only  during  and  after  the  late-winter  ice  maximum,  the  positive 

of  the  previous  summer  and  autumn  correlate  significantly  with  I  .  A 

dS 

direct  association  between  heavy  ice  and  below-normal  SSTs  is  strongly 

suggested,  especially  in  view  of  the  relatively  narrow  time  frame  (January- 

April)  during  which  the  patterns  of  southward  airflow  (P^  and  ) 

correlate  significantly  with  Bering  Sea  ice.  The  below-normal  SSTs  in 

heavy  ice  years  are  evidently  not  responses  to  the  atmospheric  flow  patterns 

that  have  contributed  to  the  heavy  ice.  The  concentration  of  significant 
HP  HP 

>  and  <IggP2  >  anto  the  winter  months  suggests  that  Bering  Sea  ice  is 


a  response  to  a  considerably  shorter  period  of  atmospheric  forcing  than  is 


19 


MP 

northern  Alaskan  ice.  The  interpretation  of  the  dual  peak  of  <1  P  >  is 

B5  X 

MP 

not  straightforward.  One  possible  explanation  is  that  positive  ,  which 

was  described  earlier  as  a  summer  mode,  is  associated  with  heavy  wintertime 

ice  via  its  strong  negative  correlation  with  summer  temperatures  in  the 

MP  P 

Bering  Sea  (positive  <P^  T^>) .  Persistence  of  these  negative  SSTs  into 

MP 

winters  of  heavy  Bering  Sea  ice  would  then  link  I  _  with  positive  P.  of 

the  preceding  summer.  This  interpretation  is  supported  by  the  tilt  of  the 
MP  P 

<P^  T^>  diagram  during  the  summer,  autumn  and  early  winter  months  (Figure  9). 

As  a  final  application  of  the  eigenvector  representation  of  the  North 

Pacific  data.  Figure  11  illustrates  the  cumulative  nature  of  the  atmospheric 

and  oceanic  for<-  .f  of  northern  Alaskan  and  Bering  Sea  ice.  The  correlations 

with  ice  extent  are  plotted  as  functions  of  the  number  of  months  (prior  to 

the  appropriate  ice  minimum/maximum)  over  which  the  atmospheric  and  oceanic 

amplitudes  have  been  averaged.  The  diagrams  therefore  show  the  dependence 

of  the  ice  indices  on  the  length  of  the  forcing  period  prior  to  31  August  in 

the  case  of  1^  and  prior  to  31  March  in  the  case  of  I  .  The  correlation 
HP 

between  I  and  P^  ,  for  example,  ranges  from  0.19  when  based  on  the  single- 

HP  HP 

month  August  value  of  P£  to  0.84  when  based  on  the  mean  value  of  P^  for 

the  9-month  period  of  December-August.  The  correlation  then  decreases  as 

the  averaging  period  is  increased  further,  i.e.,  as  months  prior  to 

HP  2 

December  are  included  in  .  The  results  imply  that  (0.84)  ~  0.71  of  the 

variance  of  the  summertime  northern  Alaskan  ice  cover  is  attributable  to 
variations  in  a  single  pattern  of  airflow  over  the  nine  months  ending  in 
August.  This  result  is  considerably  more  optimistic  than  the  forecast 
skills  of  Walsh  and  Johnson  (1979b,  Figure  18),  which  were  based  on  data 
from  single  months  and  on  a  more  general  longitudinal  formulation  of  ice 
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variability.  While  the  notion  of  the  "usefulness"  of  correlative  results  is 

subjective,  application-dependent,  and  distinct  from  statistical  significance, 

the  explanation  of  more  than  70%  of  the  variance  of  I„.  by  one  meteorological 

NA 

HP 

index,  ,  implies  considerable  potential  for  long-range  ice  forecasting  in 

the  northern  Alaskan  region.  This  potential  is  supported  by  the  physical 

HP 

plausibility  of  positive  <T  AP„  >. 

NA  2 

Figure  11  also  shows  that  the  Bering  Sea  ice  responds  to  a  relatively 

HP  MP 

shorter  period  of  atmospheric  forcing  ,  P^  )  than  does  the  northern 

Alaskan  ice,  since  the  correlations  with  I  begin  to  decay  when  the 
averaging  period  exceeds  3-5  months.  Correlations  with  the  SST  parameter 

p 

T^,  on  the  other  hand,  increase  to  0.72  at  5  months  for  northern  Alaskan 
ice  and  to  0.56  at  10  months  for  Bering  Sea  ice.  Persistence  of  the  SST 
anomalies  is  again  apparent  in  the  weaker  dependence  of  the  SST  correlations 
on  the  length  of  the  averaging  period. 

In  order  to  show  that  significant  spatial  relationships  are  not  obscured 
by  the  eigenvector  representations,  we  have  plotted  in  Figure  12  the  correla¬ 
tions  of  the  northern  Alaskan  and  Bering  Sea  ice  indices  with  actual  SLP  grid 
point  anomalies  during  the  preceding  months  of  significant  atmospheric 
forcing.  The  correlations  were  evaluated  for  all  SLP  points  on  a  5°xl0° 
latitude-longitude  grid  from  20°N  poleward.  Figure  12a  shows  that  heavy 
summertime  ice  conditions  are  preceded  by  low  pressure  to  the  east  and  high 

oressure  to  the  west  of  the  northern  Alaskan  coast.  While  I„,  correlates 

NA 

more  highly  with  the  pressures  to  the  east  rather  than  the  west,  the  highest 
7-month  correlators  are  the  pressure  anomalies  in  the  central  North  Pacific. 

A  similar  pattern  was  apparent  in  Rogers'  (1978,  Figure  2)  plot  of  ice- 
derived  summer  pressure  differences,  although  little  mention  was  made  of  the 
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Pacific  Ocean  area  in  Rogers'  work.  The  pattern  in  Figure  12a  is  suggestive 
MP  HP 

of  patterns  and  F^  in  the  North  Pacific  and  the  Arctic,  respectively. 
Interestingly,  the  highest  correlations  of  Figure  12a  are  smaller  than  those 
derived  from  the  7-month  means  of  the  eigenvector  amplitudes.  A  correspond¬ 
ing  correlation  map  for  the  3  months  preceding  the  Bering  Sea  ice  maximum 
(late  March)  is  shown  in  Figure  12b.  A  stronger-than-normal  southward 
component  of  airflow  is  indicated.  The  broad  areas  of  negative  correlations 
over  North  America  and  positive  correlations  over  the  midlatitude  North 
Atlantic  are  features  that  were  not  captured  by  the  Pacific  eigenvectors. 
However,  the  magnitudes  of  the  highest  correlations  are  similar  to  the 
eigenvector-derived  magnitudes  in  Figure  lib. 

Finally,  Figure  13  contains  maps  of  the  correlations  between  actual 
grid  point  values  of  SST  and  the  ice  indices  1^  and  IB<,.  Correlations  with 
North  Alaskan  ice  are  weakly  negative  in  the  Aleutian  region  and  more 
strongly  negative  off  the  west  coast  of  North  America.  Since  positive 

correlations  are  found  in  the  central  portion  of  the  grid,  the  overall 

p 

correlation  pattern  is  quite  similar  to  T  .  The  fact  that  1  correlates 

<v  1  JNA 

p 

more  highly  with  the  amplitude  of  T^  than  with  the  actual  SSTs  is  pertinent 

to  the  nature  of  the  I^/SST  association.  Since  most  of  the  variability 
p 

within  is  in  the  broad  region  of  positive  anomalies  in  midlatitudes 

p 

(Figure  5  ),  positive  <IN^T^>  indicative  primarily  of  the  association 
between  heavy  ice  and  positive  SST  anomalies  south  of  45°N.  This  relationship 
suggests  that  the  I^/SST  coupling  is  not  so  much  a  direct  association 
between  heavy  ice  and  low  SSTs  but  rather  an  indirect  association  involving 
the  atmosphere  (Figure  1).  The  corresponding  map  of  Bering  Sea  ice 
correlations  (Figure  13b)  contains  larger  correlations  in  the  Alaskan 
longitudes.  The  Bering  SST  anomalies 


associated  with  heavy  ice  are  negative  and  are  correlated  as  strongly  with 

p 

I  c  as  is  T, .  Again,  the  overall  pattern  of  grid  point  correlations  shows  a 
Ob  1 

P 

resemblance  to  T^.  We  conclude  that  the  eigenvector  representation  has 
captured  the  features  of  the  SST  fields  that  are  most  pertinent  to  sea  ice 
variability. 

b.  North  Atlantic  sector 

Since  the  eigenvectors  for  the  North  Atlantic  region  (Figures  6-8) 

generally  correspond  closely  with  the  Pacific  patterns  (Figures  3-5),  the 

following  discussion  will  emphasize  the  similarities  and  differences  between 

the  correlative  results  for  the  two  regions. 

Figure  14,  which  is  the  Atlantic  counterpart  of  Figure  9,  shows  the 

seasonal  and  lag  dependences  of  the  correlations  between  the  amplitudes  of 

the  North  Atlantic  SLP  and  SST  patterns.  The  diagrams  of  the  midlatitudes 

SLP/SST  correlations  show  the  same  characteristic  tilt  that  was  found  in  the 

Pacific  area.  SST  anomalies  can  again  be  interpreted  as  very  persistent 

responses  to  forcing  by  the  SLP  fields.  In  the  Atlantic  diagrams,  however, 

there  is  very  little  evidence  of  SST/SLP  associations  during  the  summer 

months.  Since  the  monthly  dependence  of  the  fraction  of  variance  described 
MA 

by  P^  does  not  show  a  non-winter  peak  of  the  type  found  with  the  correspond- 

MP 

ing  Pacific  pattern  P^  ,  there  is  indeed  no  justification  for  the  inter- 
MA 

pretation  of  P^  as  a  summer  mode.  However,  the  Atlantic  diagrams  do  show 

A  T 

that  the  second  and  third  SST  amplitudes,  T£  and  T^,  are  significantly 
MA 

correlated  with  P^  during  the  winter  half  of  the  year.  The  positive  signs 
of  these  correlations  suggest  that  the  associations  are  et  least  partially 


attributable  to  advection  by  P^  of  warm  water  into  the  western  North 
Atlantic  (T^)  and  cold  water  into  the  waters  immediately  south  of  Greenland 
and  Iceland  (T3)  • 
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Figure  14  also  shows  that  the  coupling  between  the  midlatitude  and  high- 

latitude  SLP  fields  of  the  Atlantic  region  is  similar  to  that  in  the  Pacific 

MA  HA 

sector.  Correlations  between  the  P"  and  P  are  statistically  significant, 

detectable  in  all  seasons,  and  generally  symmetric  about  0  lag.  No  systematic 

tendency  for  either  region  to  lead  the  other  is  apparent.  Similarly,  there 

are  significant  correlations  between  the  North  Atlantic  SST  fields  and  the 

high-latitude  SLP  fields.  While  the  <P^TA>  diagrams  show  a  tendency  toward 

greater  significance  at  positive  lags  of  T"A,  there  are  weak  suggestions  of  a 

feedback  in  the  opposite  direction  in  <P^AT^>  and  ‘tP^T^. 

Figure  15a  shows  the  relationships  between  summertime  ice  extent  in  the 

East  Greenland  Sea  (I™)  and  the  highest  correlators  among  the  (3-month  mean) 

EG 

patterns  of  SST  and  SLP.  The  only  significant  correlations  (r  :>  0.41)  are 

those  with  the  second,  third  and  fourth  high-latitude  circulation  patterns. 

While  even  these  correlations  are  not  large,  they  are  consistent  with 

advectional  considerations  in  the  sense  that  southward  flow  precedes  heavy 

ice.  None  of  the  midlatitude  SLP  or  SST  fields  correlate  significantly  with 

I  ,  perhaps  because  the  East  Greenland  Sea  is  the  northernmost  of  the  four 
EG 

ice  regions  in  this  study.  The  closest  grid  points  of  SST  and  midlatitude 
SLP  are  approximately  20  latitudinal  degrees  south  of  the  summertime  ice 
edge  in  the  East  Greeland  Sea. 

The  highest  correlations  with  the  Denmark  Strait  ice  index,  Ipg>  are 

show  in  Figure  15b.  The  most  striking  feature  of  Figure  15b  is  the  absence 

of  significant  correlations.  Correlations  between  and  the  SLP  amplitudes 

P^  and  P^  are  marginally  significant  (95%  level)  in  only  one  3-month 

A 

period,  December-February .  The  only  significant  SST  correlator  is  T^,  which 
correlates  negatively  with  I  and  most  strongly  in  the  months  after  the  ice 
maximum. 


As  a  counterpart  of  Figure  11,  Figure  16  shows  the  dependence  of  East 
Greenland  and  Denmark  Strait  ice  indices  on  the  length  of  the  period  over 


which  the  forcing  data  are  averaged.  The  correlations  between  I  and  the 

EG 

high-latitude  SLP  amplitudes,  and  P^\  increase  out  to  forcing  periods 

HA 

of  5-7  months.  <1  P_  >  remains  relatively  constant  until  11  months. 

EG  Z 

<IE(_,T^>  shows  a  slow  but  systematic  increase  through  11  months.  Correlations 
involving  the  Denmark  Strait  ice  index  show  little  dependence  on  the  forcing 
interval,  although  very  gradual  increases  in  the  magnitudes  of  <InQT^>  and 

i/J  J 
A 

<1  T7>  are  apparent.  <1  P  >  is  consistently  small  but  positive. 

UD  4  1/ j  / 

Plots  corresponding  to  Figure  12  were  constructed  in  order  to  illustrate 
the  relationships  between  the  North  Atlantic  ice  indices  and  the  actual  grid 
point  pressures  of  the  previous  months.  Figure  17a  shows  the  correlations 
between  t  e  East  Greenland  index  and  the  SLPs  for  the  five  months  ending  in 
July.  A  broad  area  of  negative  correlations  (pressures  below  normal  in 
heavy  ice  years)  extends  from  the  waters  east  of  the  East  Greenland  Sea  to 
approximately  170°E.  The  most  negative  values  are  near  Severnaya  Zemlya  and 
the  Laptev  Sea.  This  pattern  implies  not  only  abnormally  strong  southward 
flow  in  the  East  Greenland  passage  but  an  acceleration  of  the  entire  trans¬ 
polar  drift  stream  from  the  Chukchi  Sea  to  the  Greenland-Spitsbergen  passage. 
An  enhanced  transpolar  drift  during  the  spring  and  summer  would  be  expected 
to  contribute  to  heavy  summer  ice  conditions  in  the  East  Greenland  Sea 
(Haupt  and  Kant,  1976).  Since  Figure  17a  shows  no  high  correlations  in 
regions  other  than  the  eastern  Arctic,  East  Greenland  ice  variability  is 


evidently  decoupled  from  midlatitude  processes  to  a  greater  extent  than  is 


Figure  17b  supports  the  eigenvector-based  conclusion  that  ice  in  the 

Denmark  Strait  is  associated  only  weakly  with  the  low-level  atmospheric 

circulation  in  nearby  regions.  Marginally  significant  positive  correlations 

(above-normal  pressures  in  heavy  ice  years)  are  found  off  the  southeast 

coast  of  Greenland.  The  only  other  correlations  exceeding  0.4  in  magnitude 

MA 

are  in  the  low  latitudes,  as  suggested  by  the  negative  <1  P  >  of  Figure  15. 

Correlations  between  I  I  _  and  the  actual  grid  point  values  of  North 

ECj  ds 

Atlantic  SST  are  not  shown  because  of  the  lack  of  any  substantial  areas  of 
statistical  significance.  The  SST  grid  point  correlations  therefore  support 
Figure  15  which  showed  that  SST  patterns  do  not  correlate  well  with  subsequent 
sea  ice  in  the  North  Atlantic.  The  same  lack  of  significance  was  found  when 
ice/SST  correlations  were  recomputed  using  monthly  North  Atlantic  SST  data 
compiled  by  the  ICES  (International  Council  for  the  Exploration  of  the  Sea) 
Hydrographic  Service.  Using  the  ICES  data  set,  Schell  (1970)  described 
associations  between  SST  and  the  subsequent  atmospheric  circulation  fields 
during  selected  extreme  ice  years  in  the  North  Atlantic.  Correlations 
between  Greenland-Barents  Sea  ice  and  subsequent  SSTs  ranged  from  -0.26  to 
-0.50  for  a  set  of  49  years.  While  these  correlations  ace  com'  .-^ble  to 
those  obtained  here  (Figure  15),  Schell  did  not  present  results  for  positive 
lags  of  sea  ice.  On  the  other  hand,  Rogers  and  van  Loon  (1979)  found 
significant  correlations  between  the  atmospheric  circulation  and  subsequent 
ice  extent  in  the  Newfoundland  region,  the  Davis  Strait  and  the  Baltic  Sea. 

These  ice  fluctuations  were  associated  with  the  seesaw  of  temperature 
between  northwestern  Europe  and  northeastern  North  America.  Unfortunately, 
our  ice  regions,  the  East  Greenland  Sea  and  the  Denmark  Strait,  are  near 
the  "nodal  line"  of  the  seesaw.  The  results  from  these  previous  studies 
are  therefore  not  inconsistent  with  the  results  of  the  present  work,  although  we 


note  that  oar  conclusions  concerning  the  absence  of  significant  SST/sea  ice 
associations  at  positive  ice  lags  may  be  largely  due  to  our  particular 
choices  of  sea  ice  zones. 

Finally,  the  lack  of  significant  SLP/ice  as  well  as  SST/ice  associations 
in  the  Denmark  Strait  implies  either  that  the  ice  dynamics  in  this  region 
are  somewhat  unique,  that  I  „  is  an  inappropriate  ice  index,  or  that  the 

L)  j 

physical  signals  have  been  obscured  by  erroneous  data.  While  it  is  possible 
that  I  is  contaminated  by  some  ice  data  of  low  quality,  the  SST  consistency 
tests  described  in  Section  2  suggest  that  errors  in  the  North  Atlantic  SST 
data  are  more  likely  to  have  contributed  to  a  dilution  of  the  SST/ice 
correlations.  Nonetheless,  the  SST  and  SLP  data  show  enough  consistency 
with  each  other  (Figure  14)  that  the  small  correlations  with  IDJ,  cannot  be 
attributed  solely  to  errors  in  the  SST  and/or  SLP  data. 


5.  SUMMARY 

This  work  has  been  an  examination  of  sea  ice  fluctuations  in  the 
context  of  the  corresponding  large-scale  ice-ocean-atmosphere  systems. 

Goals  were  the  evaluation  of  the  directions  and  relative  strengths  of  the 
various  links  in  the  three-component  systems,  the  evaluation  of  the  time 
scales  of  the  various  contributions  to  interannual  ice  variability,  and  a 
comparison  of  ice-ocean-atmosphere  interactions  in  the  Pacific  and  Atlantic 
hemispheres.  While  the  results  are  limited  to  the  regions  we  have  chosen, 
the  most  pertinent  conclusions  are: 

(i)  The  dominant  patterns  of  SST  variability  are  similar  in  the 

North  Atlantic  and  the  North  Pacific.  In  both  oceans  the 
air-sea  coupling  is  such  that  the  SST  anomalies  appear  as 
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persistent  responses  to  atmospheric  forcing. 

(ii)  The  high- latitude  atmospheric  circulation  is  significantly 
correlated  with  the  circulation  fields  over  both  the 
North  Atlantic  and  the  North  Pacific.  The  midlatitude 
circulations  show  no  systematic  tendency  to  lead  or  lag  the 
arctic  circulation. 

(iii)  The  best  predictor  of  interannual  sea  ice  variability  in  the 
four  sea  ice  zones  examined  here  is  the  north-south 
component  of  the  high-latitude  circulation.  Correlations 
between  ice  extent  and  the  strength  of  this  component  over 
periods  of  several  months  are  as  large  as  0.84  in  the  mar¬ 
ginal  ice  zone  north  of  Alaska. 

(iv)  The  coupling  between  the  SST  distribution  and  subsequent  sea  ice 
extent  is  stronger  in  the  North  Pacific  than  in  the  North 
Atlantic.  The  correlative  analysis  suggests  a  direct 
associaton  between  positive  (negative)  SST  anomalies  and 
light  (heavy)  ice  in  the  Bering  Sea.  The  nature  of  the  SST/ 
ice  association  north  of  Alaska  is  most  likely  one  in  which 
the  atmosphere  serves  as  an  intermediate  link,  although  such 
a  conclusion  is  quite  tentative  in  the  absence  of  more 
detailed  information  on  ocean  current  variability  in  the 
vicinity  of  the  Bering  Strait. 

(v)  The  predictability  of  sea  ice  appears  to  be  greater  in  the  Alaskan 
sector  than  in  the  two  North  Atlantic  regions  examined  in 
this  study.  Atmospheric  forcing  over  a  period  of  up  to  9 
months  contributes  to  summer  ice  variability  along  the 
northern  Alaskan  coast,  while  atmospheric  forcing  over  3-4 


months  is  significantly  correlated  with  Bering  Sea  ice 
variability. 

(vi)  Although  sea  ice  fluctuations  in  the  East  Greenland  Sea  appear 
to  be  decoupled  from  SLP  and  SST  fluctuations  in  middle 
latitudes,  they  are  significantly  correlated  with  the 
intensity  of  the  geostrophic  airflow  parallel  to  the 
transpolar  drift  stream  during  the  previous  5-7  months. 

A  more  definitive  analysis  of  air/sea/ice  coupling  in  high  latitudes  will 
clearly  require  data  on  the  variability  of  ocean  currents  in  the  marginal 
ice  zone.  The  present  analysis  may  indeed  have  suffered  from  the  lack  of 
such  data  for  the  Bering  Strait  and  the  East  Greenland  region.  However,  our 
conclusion  that  SST/ice  associations  are  detectable  form  existing  data  in 
the  North  Pacific  but  not  in  the  North  Atlantic  underscores  the  need  for 
a  set  of  systematic  measurements  of  oceanic  variables  in  high  latitudes. 
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Correlations  between  the  sea  ice  indices  of  the  northern  Alaskan 
waters  (NA) ,  the  Bering  Sea  (BS),  the  East  Greenland  Sea  (EG), 
and  the  Denmark  Strait  (DS) .  Correlations  at  lag  0  are  computed 
from  values  for  the  same  calendar  year,  correlations  at  lag  1 
from  values  for  successive  calendar  years.  Underlined  values 
are  statistically  significant  at  95%  level. 
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FIGURE  LEGENDS 

1.  Schematic  representation  of  coupling  mechanisms  between  sea  surface 
temperature  (SST)  and  sea  ice  extent. 

2.  Departures  from  normal  ice-covered  area  (1953-1977)  in  the  northern 
Alaskan  waters  (NA)  ,  the  Bering  Sea  (BS) ,  the  East  Greenland  Sea 
(EG)  and  the  Denmark  Strait  (DS).  (See  text  for  longitudes  and 
calendar  months  corresponding  to  ice  indices.) 

3.  The  first  four  eigenvectors  of  midlatitude  North  Pacific  sea  level 
pressure.  Percentages  of  variance  described  by  eigenvectors  are 
given  in  upper  left  corners  of  diagrams. 

4.  The  first  four  eigenvectors  of  high-latitude  sea-level  pressure  in 
the  Pacific  sector. 

5.  The  first  four  eigenvectors  of  North  Pacific  sea  surface  tempera¬ 
ture.  Shaded  regions  contain  negative  anomalies  when  amplitude  of 
eigenvector  is  positive. 

6.  The  first  four  eigenvectors  of  midlatitude  North  Atlantic  sea  level 
pressure. 

7.  The  first  four  eigenvectors  of  high-latitude  sea  level  pressure  in 
the  Atlantic  sector. 

8.  The  first  four  eigenvectors  of  North  Atlantic  sea  surface  tempera- 


9.  Lagged  cross-correlations  between  amplitudes  of  dominant  eigen- 

Mt>  up 

vectors  of  mid-latitude  SLP  (P''*),  high-latitude  SLP  (P  )  and 

p 

sea-surface  temperature  (T  )  in  the  North  PaciYic.  Calendar  months 
are  identified  by  J  =  January,  F  =  February,  ...,  D  =  December. 
Contour  interval  is  0.2;  solid  contours  indicate  positive 
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Fig.  10. 


Fig.  11. 


Fig.  12. 


Fig.  13. 


Fig.  14. 
Fig.  15. 


correlations,  dashed  contours  negative  correlations.  Magnitudes 
of  correlations  exceed  0.4  in  shaded  regions. 

Cross-correlations  between  North  Pacific  seo.  ice  indices  and 
amplitudes  of  eigenvectors  of  SLP  and  SST.  Months  are  identified 
by  J  =  January,  F  =  February,  etc.  Dashed  lines  are  95%  signifi¬ 
cance  levels. 

Cross-correlation  between  North  Pacific  sea  ice  indices  and 
amplitudes  of  SLP  and  SST  averaged  over  1,  ...»  12  months  prior  to 
approximate  minimum  (northern  Alaska)  or  maximum  (Bering  Sea)  ice 
extent. 

Cross-correlations  between  (a)  1  and  grid-point  SLPs  during  7 
months  ending  31  August,  (b)  IBS  and  grid-point  SLPs  during  3 
months  ending  31  March.  Contour  interval  is  0.1.  Contours  between 
-0.4  and  +0.4  are  plotted  (dashed  lines)  only  in  region  near  North 
Pacific  ice  boundary.  H  and  L  indicate  maxima  and  minima  in 
correlation  fields. 

Cross-correlations  between  (a)  1^  and  grid  point  values  of  North 
Pacific  SST  during  5  months  ending  in  August,  (b)  TBS  and  grid 
point  values  of  North  Pacific  SST  during  5  months  ending  in  March. 
Contour  interval  is  0.2;  regions  of  negative  correlations  are 
shaded.  "+"  and  indicate  maxima  and  minima  in  correlation  fields. 

As  in  Fig.  9,  but  for  the  North  Atlantic. 

Cross-correlations  between  North  Atlantic  sea  ice  indices  and 
amplitudes  of  eigenvectors  of  SLP  and  SST.  Months  are  identified 
by  J  =  January,  F  =  February,  etc.  Dashed  lines  are  95%  signifi¬ 


cance  levels. 


A 
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Fig.  16.  Cross-correlations  between  North  Atlantic  sea  ice  indices  and 

amplitudes  of  SLP  and  SST  averaged  over  1,  ...,  12  months  prior  to 
approximate  minimum  (East  Greenland)  or  maximum  (Denmark  Strait) 
ice  extent. 

Fig.  17.  Cross-correlations  between  (a)  1^  and  grid-point  SLPs  during  5 

months  ending  31  July,  (b)  and  grid-point  SLPs  during  5  months 

ending  March  31.  Contour  interval  is  0.1.  Contours  between  -0.4 
and  +0.4  are  plotted  (dashed  lines)  only  in  regions  near  North 
Atlantic  ice  boundary.  H  and  L  indicate  maxima  and  minima  in 
correlation  fields. 


DEPARTURE  FROM  NORMAL,  105  km 


Alaskan  waters  (NA) ,  the  Bering  Sea  (BS) ,  the  East  Greenland  Sea 

(EG)  and  the  Denmark  Strait  (DS)  .  (See  text:  for  longitudes  and 
calendar  months  corresponding  to  ice  indices). 


The  first  four  eigenvectors  of  midlatitude  North  Atlantic  sea  level 


Fig.  8.  The  first  four  eigenvectors  of  North  Atlantic  sea  surface  tempera- 
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Fig.  10.  Cross-correlations  between  North  Pacific  sea  ice  indices  and  ampli¬ 
tudes  of  eigenvectors  of  SLP  and  SST.  Months  are  identified  by 
J  =  January,  F  =  February,  etc.  Dashed  lines  are  95%  significance 


levels. 
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(b)  BERING  SEA 


Fig.  11.  Cross-correlation  between  North  Pacific  sea  ice  indices  and  amplitudes 
of  SLP  and  SST  averaged  over  1,  . ...  12  months  prior  to  approximate 
minimum  (northern  Alaska)  or  maximum  (Bering  Sea)  ice  extent. 


Fig.  12.  Cross-correlations  between  (a)  1  ^  and  grid-point  SI.Fs  during  7  months 
ending  31  August,  (b)  I  and  grid-point  SLPs  during  3  months  ending 
31  March.  Contour  interval  is  0.1.  Contours  between  -0.4  and  +0.4 
are  plotted  (dashed  lines)  only  in  region  near  North  Pacific  ice 
boundary.  H  and  L  indicate  maxima  and  minima  in  correlation  fields. 


Fig.  13.  Cross-correlations  between  (a)  I  ^  and  grid  point  values  of  North 

Pacific  SST  during  5  months  ending  in  August,  (b)  1RS  and  grid  point 
values  of  North  Pacific  SST  during  3  months  ending  in  March.  Contour 
interval  is  0.2;  regions  of  negative,  correlations  are  shaded.  "+" 
and  indicate  maxima  and  minima  in  correlation  fields. 
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Fig.  16.  Cross-correlations  between  North  Atlantic  sea  ice  indices  and  ampli¬ 
tudes  of  SLP  and  SST  averaged  over  1,  ...,  12  months  prior  to 
approximate  minimum  (East  Greenland)  or  naximubi  (Denmark  Strait) 


ice  extent 


CORRELATION 


rig.  17.  Cross-correlat ions  between  (a)  I  and  grid-point  SLPs  during  5  months 

EG 

ending  31  July,  (b)  I  and  grid-point  SLPs  during  5  months  ending 
31  March.  Contour  interval  is  0.1.  Contours  between  -0.4  and 
+0.4  are  plotted  (dashed  lines)  only  in  regions  near  North  Atlantic 
ice  boundary.  H  and  L  indicate  maxima  and  minima  in  correlation 


fields. 


